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Non factorizable effects in nonleptonic B decays to charmonium 
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We discuss the validity of factorization for exclusive two-body B to charmonium transitions. In particular, we consider the role of non 
factorizable corrections in selected two-body modes. 
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Nonleptonic two-body B decays play a crucial role in the 
study of CP violation and in the determination of funda- 
mental Standard Model parameters, such as the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements or the an- 
gles of the unitarity triangle. However, it is very diffi- 
cult to deal with strong interaction effects for these purely 
hadronic processes. One of the oldest prescriptions to com- 
pute nonleptonic transitions is the factorization ansatz, of 
which various formulations have been proposed. More re- 
cent developments exploit the presence of a large param- 
eter, i.e. the b quark mass m/,. Since nib is much larger 
than the QCD scale Kqqd, it is possible to evaluate the 
relevant hadronic matrix elements as an expansion in the 
strong coupling constant a s {mi,) and in the ratio Agcz)/'«i- 
Several approaches are based on such an expansion, such 
as QCD factorization [Q] or perturbative QCD [|2). An- 
other possibility is to combine the heavy quark expansion 
with non perturbative approaches such as QCD sum rules [ 
13 or to exploit the large energy release to the final state in 
selected exclusive B decays [@J. The realm of applicability 
of each method should be assessed through the comparison 
with experimental data in each specific case. In the follow- 
ing we consider B to charmonium transitions and compare 
available experimental data to theoretical predictions based 
on factorization. In this case, such a comparison shows the 
existence of sizable violations. After a brief review of the 
factorization approach, we analyse specific decay modes 
for which the factorized amplitudes vanish; nevertheless, 
they have been observed experimentally. 



B to charmonium transitions: factoriza- 
tion versus experimental data 



Let us consider a generic nonleptonic two-body B decay: 
B — > M\M%. The effective hamiltonian describing this pro- 
cess is obtained integrating out the W field and can be writ- 
ten as: H e /f = 2,- c,(yu)(9,, where <:,(//) are short distance 
coefficients and (9, local operators. As a consequence, the 
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where /I,- represents the product of the two elements of the 
Cabibbo-Kobayashi-Maskawa matrix (CKM) involved in 
the considered transition. Long distance effects are en- 
coded in the matrix elements of the operators (9, and hence 
represent the non perturbative ingredient of the calculation. 
In order to determine such matrix elements, various ap- 
proaches have been proposed. We shall not review all of 
them here; instead, we only consider the factorization ap- 
proach (in its various formulations) with particular refer- 
ence to final products consisting of a kaon and a charmo- 
nium state. 

The effective weak hamiltonian governing the process 
B -> K~ + cc is: 
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-VtbV^CiirfOi 



i=3 



(2) 



where 0\ = (cb) V -A(sc)v-A and O2 = (sb) V - A (cc)v-A are 
current-current operators, and O3 - O\o are QCD and elec- 
troweak penguin operators [ |5) . Let us consider the mode 
B — > K~XcQ, where XcO is the ++ state of charmonium. 
Naive factorization [ 6 1 amounts to factorize the currents 
appearing in the (9, and computing the relevant amplitude 
inserting the vacuum in all possible ways. As a conse- 
quence, one has: 

A fact {B- ^ K'xco) = %V cb V* cs c(K-\(sb)v-A\B-) x 
V2 

X<X<n\(cc) V - A \Q) , (3) 

where c is a combination of Wilson coefficients. Since the 
XcO is a scalar cc particle, the latter matrix element in Q 
vanishes and therefore A f ac , - 0. Nevertheless, both Belle 
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and BaBar Collaborations have reported observation of this 
decay mode, giving: 



S(B-^^ f0 ) = (6.0!^±l.l)xlO- 4 Q, (4) 



S(B~ -> K'Xco) = (2.4 + 0.7) x 10~ 4 1 8 



(5) 



Such results clearly indicate that factorization is indeed vi- 
olated in the considered mode. Furthermore, Belle Collab- 
oration also provides the ratio: 



<B(B- -> K-J/if/) 



(0.60!o ig + 0.05 : 



: 0.08) 



(6) 



showing that the mode B 
comparable to that of B~ - 



-> K xm proceeds with a rate 
R-J/>p. 



A possible improvement of the naive factorization ap- 
proach is represented by the so called generalized factor- 
ization [|9l, in which the Wilson coefficients in the factor- 
ized amplitude are treated as free parameters. If data are 
available for some process, they can be fitted and used as 
input into other similar modes. However, it still holds that 
A(B~ — > K~Xco) = in generalized factorization. 

It is interesting to consider a decay mode allowed in this 
approximation, such as B — > K~J/if/. The factorized am- 
plitude reads: 




K 



XcO 



Figure 1. Typical rescattering diagram contributing to B 

K~Xrf>- 



(the emitted meson) is heavy, since a large overlap is ex- 
pected between Mi and the (BM\) system. An exception is 
represented by the emission of a quarkonium state, since, 
in the heavy quark limit, its transverse size becomes small. 

An analysis performed for B — > K~Xc0 shows that in- 
frared divergences in the final result do not allow to apply 
this method to such a process f 1121 . As for B — > K~J/if/, 
although the cancellation of infrared divergences has been 
proven at the leading order in the I /mi, expansion, the 
experimental data are not reproduced f 1131 . In the fol- 
lowing section, we discuss the possibility that rescattering 
diagrams, taking contribution from intermediate charmed 
mesons, could play a role in the considered modes. 



A fac1 {B- ^K-Jli/f) = 2-^V cb V; s a 2 f^x 
XFf K (Ml)(e*-q), 



(7) 



where «2 = ci + Ci/N c , is the J/ifr decay constant, q the 
kaon momentum and Ff K one of the form factors parame- 
terizing the matrix element (K~\sy M b\B~). 

Using the form factor F^ K computed in [ 10 1 and £(B~ — > 
K-J/f) = (1.01 ± 0.05) 10" 3 [ED, one obtains: a/ f = 
0.38 ± 0.05. Scanning several form factor models, the re- 
sult would vary in the range: clP = 0.2-0.4. The obtained 
effective value of «2 should be compared to the QCD cal- 
culation, which gives a^ L0 (/i = nib) = 0.163(0.126), in 
the naive dimensional regularization ('tHooft and Veltman) 
scheme [|5|. This comparison shows that non factorizable 
effects are sizable also in the case of B — > K~J/ifr. 

A QCD improved factorization approach has been pro- 
posed for B decays, exploiting the large value of mi, [ 1 1. 
The approach holds for B — > M\M2 non leptonic decays 
when M2 is light, Mi being the meson picking up the spec- 
tator quark in the decay. In this case, it has been shown that 
non factorizable corrections are dominated by hard gluon 
exchanges, while soft effects are confined to the {BM\) sys- 
tem. Naive factorization is recovered at the leading order 
in a s and KQco/mi,. The approach does not hold when M2 



3 Role of rescattering processes 

The decay B — > K~Xc0 can be obtained by rescattering 
of charmed intermediate states, as shown in fig. HI \ 1141 . 
The decay is still induced by the transition b — > sec and 
the relevant CKM structure is the same as for the direct 
transition. 

A first analysis of rescattering diagrams of the kind shown 
in fig. ffjis reported in [ 14 1 and briefly summarized below. 

The computation involves the weak matrix elements gov- 
erning the transitions B — > fl, D'*' and the strong cou- 
plings between the charmed states D { *\d ( *\ and the kaon 
and the XcO- There is experimental evidence that the cal- 
culation of the amplitude by factorization reproduces the 
main features of the B — > decay modes r 1151 . 

Therefore, neglecting the contribution of penguin opera- 
tors in d2J, we can write: 



Gf 

vr 

(DW-KV-AyO) 



{D (;y D <*)0\H w \B-) = —?=V C bV*ai{D { * )0 \(V - Af\B~) 



(8) 



with a\ = C\ + C2/N c . In the heavy quark effective theory, 
the matrix element <D ( * )0 |( V - Ay\B~) can be expressed in 
terms of a single form factor, the Isgur-Wise function £ [ 
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GD, while (D*-(p, e)|(V- A)^|0> = /d-M^, where fe is 
the D* leptonic constant and e its polarization vector; the 
analogous matrix element for D s involves the correspond- 
ing constant fo s ■ 

Other hadronic quantities involved in the calculations are 
the strong couplings of a pair of charmed mesons to the 
kaon and to the^f c o. 

The D,D^K couplings, in the soft px — * limit, can 
be related to a single effective constant g, as it turns 
out considering the effective QCD Lagrangian describ- 
ing the strong interactions between the heavy mesons and 
the octet of the light pseudoscalar mesons [ 1171 X/ = 
i g Tr[H b y,y 5 AlH a l with = ±(^-^t)^ 

iM 

and f = e J* . The matrix M contains the fields of the octet 
of the light pseudoscalar mesons. For example, the D*DK 
coupling, defined through the matrix element 



< D°(p)K-(q)\D*-(p + q, e)) >= g D .- D o K - (6 ■ q) 



(9) 



is related to the effective coupling g through the relation: 



gDTD a K- 



-2 y/m D m D - — 

JK 



(10) 



As for the coupling of the^o state to a pair of D( S ) mesons, 
defined by the matrix element 



(ID 



an estimate can be obtained considering the D( S ) matrix 
element of the scalar cc current: (Z)( S )(v')|cc|D( S )(v)), as- 
suming the dominance of the nearest resonance, i.e. the 
scalar cc state, in the (v - v') 2 -channel and using the nor- 
malization of the Isgur-Wise form factor at the zero-recoil 
point v = v'. This allows us to express gD^zWco ™ terms 
of the constant f Xc0 that parameterizes the matrix element 
(0\cc\xco(q)) = f Xc0 m Xc0 . One obtains: 



m D .,m 



gD^D, 



)L>u)Xai 



fx 



(12) 



Xm 



The method can also be applied to gjy o« ^ M . 

However, the determinations of the couplings described 
above do not account for the off-shell effect of the ex- 
changed Z)(s) and D* s) particles, the virtuality of which can 
be large. As discussed in the literature, a method to ac- 
count for such effect relies on the introduction of form fac- 
tors gi(t) = gio Fi(t), with gjo the corresponding on-shell 
couplings ( II D - A simple pole representation for Ft(f) 



is: Fj(f) = 



A z 



'DM 



Ai-t 



, consistent with QCD counting rules 



of uncertainty in our analysis. In the evaluation of the di- 
agrams, we compute at first the absorptive part, and then 
derive the real part through a dispersive representation. 

Before turning to the numerical analysis, it is worth con- 
sidering rescattering contributions of intermediate charm 
mesons to the decay mode B — > K~J/if/. The hadronic in- 
formation for determining rescattering amplitudes are the 
same as for B — > K~XcO, with the only difference in 
the strong D^D^lj/ifr couplings that can be expressed in 
terms of the parameter fj^, using the same vector meson 
dominance method applied to derive eq. dl2l . 

We have now to fix the values of the various hadronic pa- 
rameters. The Wilson coefficient a\, common to all the 
amplitudes, can be put to a\ = 1 .0 as obtained by the anal- 
ysis of exclusive B — > D s transitions. Moreover, we 
use fu, - 240 MeV, in the range quoted by the Particle 
Data Group HI II . and fjy s = fo s consistently with our ap- 
proach that exploits the large mq limit. For the Isgur-Wise 

/ 2 \ 2 

universal form factor the expression g(y) = I I is 

compatible with the current results from the semileptonic 
B D w decays. 

A discussion is needed about the D^D^K vertices. For 
the effective coupling g one can use the CLEO result g = 
0.59 + 0.01 + 0.07 obtained by the measurement of the D* 
width [ 19 1. Several estimates of g have appeared in the lit- 
erature; in particular, potential models give values close to 
one [ 20 1, while other determinations point towards lower 
values r 12 II - The discrepancy may been attributed to rel- 
ativistic effects I" 1221 . The value obtained by CLEO is in 
the upper side of the theoretical calculations f 131 1231 . We 
choose to be conservative, and vary this parameter in the 
range: 0.35 < g < 0.65 that encompasses the largest part 
of the predictions. 

In d!2t we use f XM = 5 10 ± 40 MeV obtained by a standard 
two-point QCD sum rule analysis [ 1141 . As for the cou- 
plings D^D^J/if/, expressions analogous to dl 2l i involve 
fj/t//, for which we use the experimental measurement. As- 
suming that the amplitude relative to B — > K^J/tfr devi- 
ates from the factorized result because of the contribution 
of the rescattering term: & e xp = &fact + •ftresc, one can 
constrain the values of A,- for the calculation of £(B~ — > 
K- X m)- The result is -> K'Xco) = (1.1 -3.5)x 10" 4 , 
to be compared to (1415 \ . 

Two conclusions can be drawn from the present study: i) 
rescattering amplitudes are sizeble in B — > K~xm and in 
B — > K~J/i(r; ii) they might explain the large branching 
ratio observed for B — > K'XcO- 

The calculation could be improved in several points: con- 
sidering additional intermediate states 1 , using improved 



'The contribution of orbital and radial excitations of the intermediate 
charmed mesons is expected to be suppressed by smaller values of the 



[ 1 8 1 . The parameters in the form factors represent a source leptonic constants and of the strong couplings 
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values for the input parameters or inserting new experimen- 
tal data. A more refined analysis is indeed in progress. 

It is also interesting to estimate rescattering effects in other 
channels which are also forbidden in the factorization ap- 
proximation, such as B — > K~Xc2 or B — > K h c , x& and 
h c being the 2 ++ and the l + ~ charmonium states, respec- 
tively. In particular, the latter decay mode would also be 
interesting per se, since the h c has been observed in pp an- 
nihilation, but it is not an established particle yet Nl II . A 
preliminary result gives £(B~ — > K~h c ) (1 — 3) X 10~ 4 [ 
24 1 to be compared to the prediction for the inclusive mode 
S(B~ -> Xh c ) (0.13 - 0.14)% 

It has been suggested that a possible decay chain to ob- 
serve this mode could be: B — > K~h c , h c — > yr\ c with 
T] c — > KKn or ij c — > t/ww, and S(/i c — > 77^7) - 0.50 + 0.11 
r 1261 . Using such a prediction, our result would give: 
S(B- -> K h c -> Kr] c y -> K(Kkn)y) * (2.5-7.5)x 10~ 6 , 
suggesting that this mode could be accessible at the cur- 
rently operating experimental facilities. 

4 Conclusions 

Present experimental results on two-body B to charmo- 
nium transitions show large non factorizable contributions. 
We suggested that such contributions can be interpreted by 
rescattering of intermediate charmed resonances. The nu- 
merical evaluation of such effects for B — > K~Xc0 shows 
some agreement with experimental data. Contributions of 
similar sizes are expected for B — > K ~h c . Predictions for 
other decay modes could confirm this picture. 

Partial support from the EC Contract No. HPRN-CT-2002- 
00311 (EURIDICE) is acknowledged. P. Colangelo and 
T.N. Pham are thanked for fruitful collaboration on the top- 
ics discussed above. 
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